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The spectra of lepton pairs from correlated open charm and bottom decays in ultrarelativistic 
heavy-ion collisions are calculated. Our approach includes energy loss effects of the fast heavy 
CN ■ quarks in deconfined matter which are determined by temperature and density of the expanding 
fvj . parton medium. We find a noticeable suppression of the initial transverse momentum spectrum 
of heavy quarks due to the energy loss within the central rapidity covered by the PHENIX 
Q^ . detector system. As a result the dilepton yields from correlated decays of heavy quarks are in 
the same order of magnitude as the Drell-Yan signal. 

2 1 Introduction 

The production of high invariant mass dileptons in relativistic heavy-ion collisions has 
rS I attracted a great deal of interest in the recent years from both the experimental and 
5^ I theoretical sides. Dilepton measurements are planned with the PHENIX and ALICE 
detector facilities at the relativistic heavy-ion collider RHIC in Brookhaven and at the 
large hadron collider LHC in CERN, respectively. In despite of many theoretical attempts 
the nature of the dominating source of dileptons is still unclear. Among the candidates 
that can give strong contributions in the large invariant mass continuum region near and 
beyond the J/ip one can consider the following processes in which dileptons are generated: 
(i) hard initial quark - anti-quark collisions give rise to the Drell-Yan yield, (ii) semi- 
leptonic decays of heavy quarks like the charm and bottom, which are also produced 
dominantly in hard initial parton interactions, and (iii) the thermal dilepton radiation 
resulting in interactions between secondary partons or even in a locally thermalized quark- 
gluon plasma (QGP). The yield from (iii) is expected to serve as a direct probe from the 
QGPjj. 
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It is commonly believed that, due to the large invariant mass M, at least the Drell- 
Yan process and the initial heavy quark production in pp and pA collisions are under 
reliable theoretical control by means of perturbative QCD and, therefore, can provide a 
reference for the thermal dilepton signal from deconfined matter. At the same time all up 
to date calculations of high-M thermal dilepton rates from early parton matter contain 
uncertainties due to the lack of knowledge of precise initial condition parameters such as 
the energy density and quark-gluon phase space saturation. Recent estimates 0] within 
the mini-jet mechanism of quark-gluon matter formation Q, with initial conditions similar 
to the self-screened parton cascade model 0, point to a strong competition between the 
Drell-Yan yield and the thermal signal in a wide range of beam energies. 

With respect to the wanted dilepton signal from the QGP the dileptons from charm 
and bottom decays play a role of substantial background |0. The magnitude and the 
characteristic properties of such background is presently matter of debate. As shown by 
Shuryak the energy loss effects of heavy quarks propagating through a medium can 
drastically reduce the dilepton yield from charm and bottom decays in the region of large 
invariant dilepton masses. Further calculations 0, including the PHENIX acceptance at 
RHIC, also predict a very strong suppression of dileptons from heavy quark decays, which 
then is partially below the Drell-Yan yield. The studies in refs. [^ ^ are based on the 
approximation of a constant rate of energy loss —dE/dx = 1 or 2 GeV/fm. However, 
the value of dE/dx is expected to depend on the mean free path of heavy quarks in 
the medium which is obviously governed a time dependent density. As shown below, in 
particular within the gluon radiation mechanism of energy loss 0, dE/dx can depend on 
widely varying characteristic quantities of the parton medium. 

In the present contribution we study the dileptons resulting from decays of correlated 
open charm and bottom produced initially in hard collisions and undergoing an energy 
loss which is determined by the temperature and the density evolution of the deconfined 
medium. 

2 Initial hard production processes 

In the framework of perturbative QCD the number of the heavy quark - anti-quark (QQ) 
pairs, produced initially with transverse momenta p±i = —p±2 = P± at rapidities j/i and 
1/2, in central AA collisions can be calculated from 

d^QQ = Taa(O) ICq H{yi, y2,p±) dp\ dyi dy2, (1) 

where H{yi,y2,p±) is the standard combination of structure functions and elementary 
cross sections (see for details [01 |6|) 



H{yi,y2,p±) = xiX2<fg{xi,Q'^)fg{x2,Q'^'^ — - 
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+ j:[f,ixi,Q')Mx2,Q') + f,ix2,Q')Mx„Q')]^\, (2) 



where fi{x,Q) with i = g,q,q are the parton structure functions, a;i,2 



mj_ (exp{±|/i} + exp{±y2}) /y/s and mj_ = Jpj_ + itiq. As heavy quark masses we take 
rric =1.5 GeV and mb = 4.5 GeV. We employ throughout the present paper the HERA 
supported structure function set MRS D'- from the PDFLIB in CERN. The overlap func- 
tion for central colhsions is T4a(0) = A'^/{TrR\) with Ra = I.IA^^^ fm and A = 200 in this 
paper. We do not include shadowing effects of nuclear parton distributions in the present 
paper since for heavy quark production they are expected to be not very important and 
can be considered separately. Our calculation procedure is based on the lowest-order QCD 
cross sections da^g/dt for the subprocesses gg -^ QQ and qq —>■ QQ with the simulation 
of higher order corrections by the corresponding /Cq factor. Such a procedure reproduces 
within the needed accuracy the more involved next-to-leading order calculations ^ for 
the heavy quark pair distributions with respect to their invariant mass, total pair rapidity 
and gap rapidity. We find the scale Q^ = 4mQ and /Cg = 2 as most appropriate. 

For the Drell-Yan production process of leptons at rapidities yi and 1/2 and transverse 
momenta p^i = "~P±2 = Pi. we have 



dN^j = TAA{0)ICDYL{yi,y2,P±)dpj_dyidy2, (3) 

L{yi,y2,pi_) = ^xiX2 fgixi,Q^) fg{x2,Q'^) + fgix2,Q'^) fgixi,Q'^) — |-, (4) 
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dt 



with da^^/di = ^ch{yi - 1/2)/ (2pl[l + ch(yi - y2)f), Xi,2 = p± (exp{±|/i} + exp{±|/2}) 
/y/s, a = 1/137, Q^ = X1X2S = M^ and ICdy = l-l- 

3 Propagation of the heavy quarks through the par- 
ton medium and energy loss 

Since the energy loss of heavy quarks depends strongly on the properties of produced 
matter, which is propagated through, one has to specify the initial conditions and the 
space-time evolution of such a medium. To do this in a coherent manner we employ for the 
mini-jet production (which dominates the parton matter formation) the same lowest-order 
approximation and parton structure functions as for the initial heavy quark production. 
Adding a suitably parametrized soft component ||^ we get the initial temperature T, = 
550 (1000) MeV, gluon fugacity Xf = 0.5 and light quark fugacity A^ = Xf/5 of the mini- 
jet plasma formed at initial time Tj = 0.2 fm/c at RHIC (LHC). The space-time evolution 
of the produced parton matter after Tj is governed by the longitudinal scaling-invariant 



expansion accompanied by quark and gluon chemical equilibration processes |T0|, ^^ . We 
take for definiteness full saturation, i.e. A^ = A'' = 1 at deconfinement temperature 
Tc = 170 MeV. The actual time evolution is determined by e oc r~^/^ for energy density 
and A^'^ocr^ Q. 

^ Due to gluon multiplication this seems to be justified, while the chemical equilibration of the quark 
component might be somewhat overestimated by this ansatz. 



To model the energy loss effects of heavy quarks in expanding matter we assume 
as usual in the Bjorken scaling scenario that a heavy quark produced initially at given 
rapidity follows the longitudinal flow with the same rapidity. Therefore, with respect to 
the fluid's local rest frame the heavy quark has essentially only a transverse momentum 
P±(t) depending on the proper local time r in accordance with the energy loss in transverse 
direction. 

To calculate the evolution of the transverse momentum P±{t) of quarks propagating 
the distance r_|_(r) in transverse direction we adopt the results of Baier et al. [^ for the 
energy loss of a fast parton in a hot QCD medium 



dp± 4a,fcc,2 , 2U/4 1 [VPj- + "^Qi /.N 

^ = -371^''" + "^^ ^"1, Lfcg I' (') 

dr± pj_ 

dr Jpl + ml' 



(6) 



where the parameter k^ = 2mff^(T) = |(Ag + jXq)TcasT^ [T^, |n| is used as an average 
of the momentum transferred in heavy quark-parton scatterings, and the strong coupling 
is described by a^ = 0.3. For the mean free path of gluons we take L ^ {.^gg'<^g)~^ with 
integrated cross section agg ~ (^1/k1 and gluon density Ug ~ Ac,T^, resulting in L~^ = 
2.2 a^T |jTO|. The derivation of eqs. (||,|) relies on the assumption that the temperature 



changes smoothly enough in space-time so that the mean free path of heavy quarks is less 
than the characteristic scale of the temperature gradient. Such an approximation can be 



improved by considering finite size effects |T2[. The r.h.s. of eq. (^) is actually the rate of 
energy loss dE/dx which is essentially oc {p^ + ?71q)^/^. Due to the rapid expansion of the 
matter this value depends strongly on the time. For instance, at RHIC conditions and 
initial momentum of a heavy quark of 5 GeV the stopping power —dE/dx drops from 2.8 
GeV/fm down to 0.5 GeV/fm during the expansion of deconfined matter. 

To get realistic spectra of charm and bottom quarks after energy loss we use a Monte 
Carlo simulation with a uniform distribution of the random initial position and random 
orientation of QQ pairs in the transverse plane. We integrate then eqs. (|^,§) together 
with the above described time evolution of T{t) and Xg^g^r). Eq. (|^) is used to check 
whether the considered heavy quark propagates still within deconfined matter; if it leaves 
the deconfined medium it does not longer experience the energy loss according to eq. (^). 
(The subsequent energy loss in a possible mixed and hadron phase might be incorporated 
along the approach of |]T3| •) A considerable part of heavy quarks can not escape the parton 
system thereby undergoing thermalization. Following [§ we consider a heavy quark to be 
thermalized if its transverse mass m± during energy loss becomes less than the averaged 
transverse mass of thermalized light partons at given temperature. This part of heavy 
quarks we assume to be distributed as dNQ/dmj_ oc m^ exp{— mj_/Tcfr} with an effective 
temperature T^s = 150 MeV. 

The results of our calculation for transverse momentum spectra of charm and bottom 
quarks at RHIC are shown in fig. 1. We obtain a noticeable suppression of the initial 



transverse momentum spectrum at large transverse momenta due to the energy loss both 
for charm and bottom. The bump below 3 GeV is caused by the thermal redistribution. 



In comparison with LHC conditions |]I4[ the suppression at RHIC is found to be less pro- 
nounced due to the shorter life time of deconfined matter and its lower initial temperature 
and density. It should be stressed that such a change of the p± spectrum of heavy quarks 
in AA collisions, in comparison with pp or pA collisions, is a measurable effect which for 
itself can help to identify the creation of a hot and dense parton gas [1^ . 



4 Dilepton spectra from heavy quark decays 

Basing on the analysis of ref. P we employ a 6 like fragmentation function for heavy 
quark fragmentation into D and B mesons. The correlated lepton pairs from DD and BB 
decays have been obtained from a suitable Monte Carlo code which employs the single 
electron and positron distributions as delivered by JETSET 7.4. The average branching 
ratio of D ^ l+X is taken to be 12%. We consider here only the so-called primary leptons 
PI which are directly produced in the bottom decay b —>■ I X with the average branching 
ratio 10%. The secondary leptons /' produced in further branchings like b ^ cX -^ I'Y 
have much smaller energy and do not contribute to the large invariant mass region. Such 
simplification does not affect significantly the dielectron yield at M > 2 GeV as we 
have checked by comparison with ref. [^]. To get the dielectron spectra from charm and 
bottom decays at RHIC we take into account the acceptance of the PHENIX detector: the 
electron pseudo-rapidity is |?7e| < 0.35 with an additional cut on the electron transverse 
momentum p']_> 1 GeV to reduce the lepton background from light hadron decays. Our 
results of heavy quark decays are depicted in fig. 2 for charm (a) and bottom (b) for 
M > 2 GeV. Clearly seen is that both the charm and bottom signals are suppressed 
so that they compete with the Drell-Yan yield. Changing the kinematical gates (i.e., 
narrowing rapidity or enlarging the minimum p±) does not help to suppress the charm 
and bottom signals below Drell-Yan. 

To demonstrate the general tendency of the relative contribution from the Drell-Yan 
process and heavy quark decays with the change of the collider energy we plot in fig. 3 
the results of our calculations for LHC initial conditions including the ALICE detector 
acceptance for the electrons, i.e. \rie\ < 0.9. As above we here also employ the gate 
p5_ > 1 GeV. One observes that the strong energy loss at LHC is not enough to suppress 
sufficiently the signal from bottom decays, while the charm yield is clearly below the 
Drell-Yan yield. For more details see |jl^ . 



5 Summary 

In summary, taking into account the energy loss effects of heavy quarks we calculate the 
spectra of electron pairs from correlated open charm and bottom decays. We obtain a 
measurable suppression of large transverse momenta of heavy quarks due to the energy loss 



in the expanding deconfined matter. Our calculations predict competing contributions of 
charm and bottom decays and Drell- Yield in the high invariant mass region of electron 
pairs which are planned to be measured at central rapidity in the PHENIX detector. 
Further studies, including a different fragmentation scheme and other kinematical cuts, 
are in progress. 

Acknowledgement: Stimulating discussions with participants of the RHIC'97 work- 
shop are gratefully acknowledged. 
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Figure 1: The p± spectrum of heavy quarks ((a): charm, (b): bottom) in central AA 
colhsions at RHIC. The dashed curves depict the initial production without energy loss, 
while the solid lines show our result with energy loss and Tcs = 150 MeV for thermalized 
heavy quarks. Fat curves are for full phase space, while the light curves display these heavy 
quark pairs which are detected in PHENIX through their primary decay dielectrons. 
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Figure 2: The invariant mass spectrum of dielectrons from correlated charm (a) and 
bottom (b) decays as well as Drell-Yan electron pairs filtered throughout the PHENIX 
acceptance. 
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Figure 3: The invariant mass spectrum of dielectrons from correlated charm (short- 
dashed hne) and bottom (sohd hne) decays after enegy loss as well as Drell-Yan electron 
pairs (long-dashed line) filtered throughout the ALICE acceptance as described in text. 



